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I. INTRODUCTION 

In communication theory, it is generally believed that mem- 
ory effects improve the information transfer capabilities of a 
communication line Memory effects can be introduced 
both by the presence of correlations in the noise affecting 
different channel uses (inputs), and by interference among 
theme. 

In the quantum communication scenario the problem of de- 
termining the optimal ensemble of input states, depending on 
the memory, naturally arises. An optimal ensemble is the 
most robust under the action of the noisy channel, leading to 
the highest transmission rates. In particular, for the classi- 
cal information transmission through quantum memory chan- 
nels, the possibility of discriminating between the optimality 
of separable or entangled input states would be useful. 

Then, by the transitional behavior it is intended the possi- 
bility to single out two "phases" in the channel capacity. One 
for which the optimal input states are entangled among dif- 
ferent channel uses, and the other for which the optimal in- 
put states are separable. To model the memory effects it is 
customary to introduce a memory parameter, which quanti- 
fies the amount of correlations in the quantum channel. The 
transition between the two "phases" may happen at a finite 
value of the memory parameter, implying that separable in- 
put states are optimal in the presence of small correlations, or 
at zero value of the memory parameter, implying that sepa- 
rable states are optimal only in the memoryless limit. Sev- 
eral models showing such an effect have been proposed for 
discrete quantum memory channels [2] as well as for contin- 
uous ones [3]. However, the majority of these works restrict 
their analysis to few channel uses (thus not relying on capacity 
arguments) and above all do not provide general criterion to 
characterize the transitional behavior. Here we present neces- 
sary conditions to have the transitional behavior for a class of 
quantum Gaussian memory channels. Actually we show that 
such behavior is intimately related to the symmetry properties 
of the channels' action. This is possible thanks to the intro- 
duction of a technique which allows us to unravel the memory 
effects. 

The article is organized as follows. In Sec. HI] we review 
the basic tools for working out quantum Gaussian memory 
channels, and to the problem of evaluating the channel ca- 
pacity. In Sec. Hn] we consider the problem of computing 
the Holevo function for a Gaussian memory channel, and in- 



troduce a class of Gaussian memory channels for which this 
problem can be traced back to the memoryless setting. In Sec. 
llVlwe enunciate a criterion for the transitional behavior and 
we relate the optimality of entangled inputs to the symme- 
try properties of the channels' action. In Sec. IVl illuminating 
examples are presented. Finally, Sec.|Vl]is for concluding re- 
marks. 



II. GAUSSIAN STATES AND GAUSSIAN CHANNELS 

Gaussian quantum states and Gaussian quantum channels 
are defined in the context of continuous variable quantum 
systems. Here we consider the case of a continuous vari- 
able quantum system consisting of n identical quantum har- 
monic oscillators (for a complete presentation of the subject 
see for instance |4]). The n quantum harmonic oscillators 
are associated with a set of 2n canonical operators qi,pi, 
<l2,P2, ■ ■ ■ qn,Pn, Satisfying canonical commutation relations 
[qh,Pk] = i^hk (here and in what follows we assume h = 1). 

The state of the n modes can be described, using the for- 
malism of density operator, as a certain density /5„ defined in 
the n-mode Fock space. However, we find it more convenient 
to work in the Wigner function representation. Let us recall 
that the Wigner function associated to a density p„ is defined 
as 



W"n(q,p)= / d"y(q-y|/5„|q + y)e 



2i'7ry -p 



(1) 



where we have introduced the numerical vectors q := 
{qi, . ..qn), p [pi, . . .p„), y := (yi, ...yn) e K", and 
we have denoted 



|q±y) ■.= ®l^^\qk±yk) 



(2) 



the joint (generalized) eigenstates of the 'position' operators 
{'7fe}fc=i,...n, i-e. gfc|q ± y) = {qk ± ?;fe)|q± y>, for all k = 
1, . . .n. 

By definition, Gaussian states are those described by a 
Gaussian Wigner function. In the ?i-mode scenario, the 
Wigner function of a Gaussian state is a multivariate Gaus- 
sian function: 

Ty„ X = ^ ^, (3) 

^ ' (2^)"v/dit(T0 
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where we have introduced the numerical vector x := 
iqi,Pi,q2,P2, ■ ■ ■qn,PnV e R^"- The Wigner function of 
an n-mode Gaussian state is hence completely described by 
the vector of the 2n first moments 

m= (x) = y"xVK„(x)d"x, (4) 

and by the 2n x 2n covariance matrix (CM) 

V = ((x-m)(x-m)"^) = J (x-m)(x-m)"^T4^„(x)d"x. 

(5) 

Finally, let us notice that certain conditions have to be im- 
posed on the CM, to ensure that the Wigner function describes 
a bona fide quantum state. Indeed, the Heisenberg principle 
imposes the condition f3] 

V -in>0, (6) 

where 

fe=i ^ ^ 

is the 2n x 2n matrix representing the n-mode symplectic 
form. 

We will also consider an equivalent representation defined 
by a different ordering of the canonical variables, expressed 
by the numerical vector X = {qi,q2, . . .qn,Pi,P2, ■ ■ ■Pn)~^ , 
with the corresponding vector of first moments m and the CM 

1/= ((x-m)(x-m)"r). (8) 

In this representation, the symplectic form is represented by 
the matrix 

where O denotes the null matrix, and I„ the identity matrix of 
size 71. 

A Gaussian quantum channel acting on n bosonic modes (in 
short, an ?i-mode Gaussian channel) is by definition a channel 
mapping Gaussian states into Gaussian states. As a conse- 
quence, its action on Gaussian states is completely character- 
ized by the rule of transformations of the vector of first mo- 
ments and of the CM. One can show (see, e.g., |i^) that a 
Gaussian channel transforms the pair (m, V) (vector of first 
moment, CM) as follows: 

{m,V) ^ iXm + d,XVX'^ + Y). (10) 

Where d e R-^" is a displacement vector, and X, Y are two 
2n X 2n matrices. In order to represent a bona fide quantum 
channel, these matrices have to obey the inequalities 

Y + ixnx'^ -in>o, Y>o, (i i) 

In conclusion, a Gaussian channel is characterized by the triad 
(d, X, Y), satisfying Eq. 



Given a pair of Gaussian channels: $ with associated triad 
(d, X, Y), and with triad (d', X' , Y'), the composition of 
the two channels o $ is associated to the triad {X'd + 
d',X'X,X'YX'~^ + Y'). 

In the family of Gaussian channels, a special sub-family is 
those of unitary Gaussian channels. Unitary Gaussian chan- 
nels are characterized by the conditions Y — 0, and 

Xnx'^ = n. (12) 

The last equation characterizes linear symplectic transforma- 
tions, where the matrix X is symplectic. In conclusion, a 
Gaussian unitary transformation is characterized by the triad 
(d, X, Y), where d is generic, X is symplectic, and Y — 0. 
In the following we consider the subgroup of Gaussian uni- 
tary transformations preserving the total number of excita- 
tions, such maps are represented by matrices X which are 
both symplectic and orthogonal. Using the representation ([8]l, 
it is possible to show (see, e.g., |4] and the references therein) 
that those matrices are of the form 




where A, B are real matrices of size n such that the matrix 
A + iB is unitary. 

A. Gaussian memory channels, normal forms 

For any integer n, n uses of a Gaussian channel transform 
n input bosonic modes into n output bosonic modes. The ac- 
tion of a Gaussian quantum channel is hence described by 
a sequence of Gaussian channels $„, acting on n bosonic 
modes, which is in turn associated to the sequence of triads 

(dn 1 X^ , Yji 'j. 

A very special case is that of the memoryless channels, for 
which <&„ = $f " is the direct product of n identical one- 
mode Gaussian channels. Each of these identical one-mode 
Gaussian channels is characterized by a triad (di,Xi, Yi). 
Hence, a memoryless channel is characterized by a sequence 
(d„,X„,y„) = (eLidi,eLi^i,©Lii"i)- Notice 
that di — {dq^dpY, and we have denoted di := 

{dq,dp, dq, dp, . . . dq, dpY e K^". Also notice that Xi, Yi 
are 2 X 2 matrices, and the direct sums Xi, 
are 2n x 2n block-diagonal matrices. 

Let us now consider the case of a quantum channel with 
memory. We call a channel with memory, or simply a memory 
channel, any channel which is not memoryless. Making no as- 
sumption on additional structures that might be present (e.g. 
causality, invariance under time translations), we can only say 
that the associated sequence of n-mode Gaussian channels 
satisfies $„ ^ (j)®"^ and the associated sequence of triads 

is (d„, y„) ^ (e;:^, di, e;'^, ^i, ©Li ^i)- 

The problem of finding normal forms for n-mode Gaus- 
sian channels was considered in , the case of one-mode 
channel was considered in |9]. Normal forms are equivalence 
classes of n-mode Gaussian channels, up to (Gaussian) uni- 
tary equivalence. Hence, given a pair of n-mode Gaussian 
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channels $„, they are equivalent if there are Gaussian 
unitary transformations £„, I?„ such that — I?„ o $„ o f „ 
(o denotes the composition of channels). As the chosen nota- 
tion suggests, we want to look at £„, I?ri respectively as uni- 
tary encoding and decoding transformations, £„ anticipating 
and 2?„ following the action of the quantum channel. We will 
hence write (d„, Xn, Yn) ~ (d'j, X^^,YJ^) if the two Gaussian 
quantum channels are equivalent in the sense declared above. 

The first thing to be noticed is that {dn,Xn,Yn) ~ 
{0, Xn,Yn). This is a well known result, a consequence of 
the fact that the displacement vector can always be eliminated 
by applying a proper n-mode displacement operator (which is 
a Gaussian unitary transformation) at the encoding, or decod- 
ing, stage. For this reason, in what follows we only consider 
71-mode Gaussian channels of the form {0, Xn,Yn). Let us 
consider a quantum channel $„ with the triad (0, X„, an 
encoding £„ with triad (0, En, 0), and decoding 2?„ with triad 
(0, Dn,0)- The application of the encoding and decoding uni- 
taries leads to the dressed channel $'j associated to the triad 

{0,Dr,X„En,DnYnl^). 

As was shown in Q\, an n-mode Gaussian channel is al- 
ways unitary equivalent to a n-mode channel in a normal 
form, i.e. (0, X„, y„) ~ (0, X'^, F^), where 



X' — 



(h) 

2 



L/i=i 



e 



fc=2p+l 



(14) 



is the direct sum of p two-mode (4 x 4) matrices and 

n — 2p one-mode (2 x 2) matrices X^ . In other words, by 
applying suitable encoding and decoding unitaries, the matrix 
Xn is reduced to the direct sum of two-mode and one-mode 
terms. However, it is important to notice that the matrix 
cannot be in general jointly reduced to the same form. 



B. Classical capacity of Gaussian channels 

A quantum channel can be used to transmit classical infor- 
mation by encoding a classical stochastic continuous variable 
Z, distributed according to a probability density distribution 
pz, into a set of quantum states pz- 

In the case of a memoryless quantum channel $„ = cj)®"^ 
the maximum rate at which classical information can be reli- 
ably sent through the quantum channel is given by the regu- 
larized limit 1,10.1 

C = lim -X i^T) , (15) 
where the Holevo function x evaluated on n channel uses is 



max <, S 

{pz,pz} 



dZpzpz 



-J [$f " (pz)] I , (16) 



where S denotes the von Neumann entropy. If the von Neu- 
mann entropy is expressed in qubits, then the classical capac- 
ity of the quantum channel in measured in bits per channel 



use. The computation of the memoryless channel capacity is 
based on the optimization over all input ensembles, includ- 
ing those made of states which are entangled among different 
channel uses. On the other hand, if the input states are re- 
stricted to ensembles of separable states, one obtains the so- 
called one-shot capacity 



Ci^xi'^i) 



(17) 



Clearly the one-shot capacity is a lower bound on the mem- 
oryless channel capacity. If the two quantities coincide, the 
Holevo function is said to be additive. Additivity of the 
Holevo function dramatically simplifies the problem of eval- 
uating the memoryless channel capacity. Even though the 
Holevo function has been shown to be additive for several 
relevant channels, e.g. the lossy channel in the framework of 
Gaussian channels ifllll . this property does not hold in general 
1 12]. Non-additivity of the Holevo function implies that the 
optimal ensembles of input states, i.e. the most robust to noise, 
are entangled among different channel uses: a phenomenon 
which has no counterpart in the classical theory of informa- 
tion. 

Moving to the case of quantum channels with memory, 
characterized by the inequality $„ ^ (j,®"^ Qjjg could be 
tempted to generalize the formula in Eq. ( fT5] l and write 



C: 



lim -xi^n) , 



(18) 



with 



X($„)= max is 

{pz,Pz} 



<!>„, I / dZpzpz 

dZpzSi^nipz)] 



(19) 



Indeed, it is possible to show II 1311 that the right hand site 
of Eq. ( fTSl l is in general only an upper bound for the classi- 
cal capacity of the memory channel. On the other hand, it 
has been proven that this quantity coincides with the mem- 
ory channel capacity for the class of so-called forgetful chan- 
nels ifTill . Those channels have the property that correlations 
among channels uses decay exponentially. Moreover, they ex- 
hibit a causal structure and are invariant under time transla- 
tion. 

The problem of computing the classical capacity is ex- 
tremely hard. Indeed, if one cannot rely on the additiv- 
ity property, one has to evaluate the regularized limit of the 
Holevo information, whose complexity increases exponen- 
tially in n. Moreover, for the case of bosonic channels, the 
relevant Hilbert space is infinite dimensional even for a sin- 
gle channel use, i.e. n ~ 1. Clearly, an infinite dimensional 
Hilbert space could carry an infinite amount of classical bits. 
Hence, to avoid unphysical results, one is led to introduce a 
physically motivated constraint, and to exploit the concept of 
constrained capacity. A typical choice in the framework of 
bosonic Gaussian channels is to impose a constraint on the 
maximal input energy per channel use, hence one introduces 
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the constrained Holevo function 



X" ($») = max 

{pz.Pz} 



E 

fe=l 



s 



tr 



dZpzpz 
dZpzS [$„ (p^)] 



91 



(20) 



where we have assumed unit frequency for the bosonic os- 
cillators, and N represents the maximum number of excita- 
tions per mode on average. One can conjecture that the max- 
imum is reached in correspondence with a Gaussian ensem- 
bles. Indeed, the optimality of the Gaussian ensemble has 
been proven for the lossy channel lITTIl and conjectured for 
other families of bosonic Gaussian channels ifTsll . 

Here, we estimate the Holevo function when restricted to 
Gaussian encoding 1 16]. We consider Gaussian encoding de- 
fined as follows. For n uses of the quantum channel, we fix a 
reference n-mode Gaussian state, with zero mean, described 
by the Wigner function 



exp[-lxTx/-ix] 
(27r)"v/det(V^„) 



(21) 



A classical variable m G K^" is hence encoded by applying a 
displacement operation on the reference state, thus obtaining 



exp[-l(x-m)Tl^-i(x-m)] 
(2^)"v/det(T/i„) 



(22) 



We assume the stochastic variable m to be itself distributed 
according to the Gaussian probability density distribution with 
zero mean: 



exp [— im^Vj. ^m] 
(2^)»v/drt(^ 



(23) 



By linearity of the relation (HJ, the corresponding ensemble 
state 



/ 



is itself described by a Gaussian Wigner function, i.e. 



exp[-^xT(T/i„ + -l/)-ix] 
(2^)"Vdet(Fin + yc) 



(24) 



(25) 



The restriction to Gaussian states, which are mapped into 
Gaussian states by Gaussian channels, dramatically simpli- 
fies the problem, since the complexity of Gaussian states is 
polynomial in the number of modes n. Moreover, the von 
Neumann entropy can be calculated in terms of the symplec- 
tic eigenvalues of the CM. 

The symplectic eigenvalues of an n-mode CM V are de- 
fined as follows. Notice that the matrix Vfl, where fl is the 



symplectic form introduced in Eq. (O, has 2n purely imag- 
inary eigenvalues {±ii'k}k=i....n, where the n real numbers 
{i'k}k=i,...n are the symplectic eigenvalues of the CM V. 
From the uncertainty relations, expressed by Eq. it fol- 
lows that the symplectic eigenvalues satisfy the inequalities 
i^k > 1/2, which are saturated by pure Gaussian states 101 ■ 
The von Neumann entropy of an n-mode Gaussian state p, 
characterized by a CM V, is given by the formula; 



S[p] = ^9\'^k 

k=l 



(26) 



where we have introduced the function g defined by 

g{x) = (x + 1) log2 {x + l)-x log2 (x) . (27) 

The von Neumann entropy of a Gaussian state p is determined 
by its CM V, it is hence convenient to define a function S of 
the CM such that 



nv] = s{p) . 



(28) 



Finally, we notice that the input energy constraint in Eq. 
(|20] i can be written in terms of the CM as follows: 

2n 2 



III. THE HOLEVO FUNCTION FOR GAUSSIAN 
ENSEMBLES 



For n uses of the quantum channel, the constrained Holevo 
function, when restricted over Gaussian ensembles, reads 

Xci^n) = max - {S [X„(M„ + Vc)Xj + r„] 

Vin.Va n 



tr(T4n + Vc) ^ 1 
2n - 2 



(30) 



where the optimization is over the CMs Vin, K satisfying the 
energy constrains. 

In the case of a memoryless Gaussian channel, by restrict- 
ing on Gaussian input states which are separable among dif- 
ferent channel uses, we get to the one-shot capacity: 

X^($i) = max {S [XiiVin + Vc)Xj + Y^] 



- S [XiVinXj + Yi] 

2 - 2 



(31) 



In general, for a Gaussian memory channel, characterized 
by the sequence of triads (d„,X„,y„), the optimization of 
the n-use Holevo function in Eq. ( [30l l cannot be reduced to the 
one-use case in Eq. (|3T| i. That is a consequence of the the nor- 
mal form for an n-mode Gaussian channel in Eq. ( fT4] i. which 
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in general is not in the form the product of n independent one- 
mode Gaussian channels. However, we can still identify a 
class of Gaussian memory channels such that, for any n, there 
are Gaussian unitary encoding and decoding transformations 

£n = (0, En, 0), Vn = (0, Dn,0), such that 



o $„ o = (g) 



(k) 



(32) 



k=l 



where {^^''^Ife^i,. „ is a collection of n (not necessarily iden- 
tical) one-mode Gaussian channels. In other words, the mem- 
ory channels belonging to that class factorize in terms of the 
collective input and output variables defined by the encoding 
and decoding transformations. Hence we introduce the fol- 
lowing 

Definition 1 (memory unraveling) A bosonic Gaussian 
memory channel, characterized by a sequence (0,X„,y„), 
can be unraveled if there is a sequence of encoding Gaussian 
unitaries {Q,EmQ), and a sequence of decoding Gaussian 
unitaries (0, D„, 0), such that, for any n. 



D Y 



k=l 
n 



(33) 
(34) 



k=l 



Moreover, we require that the encoding unitary preserves the 
form of the energy constraint. 



Using this ansatz, the calculation of the Holevo function for n 
uses of the memory channel reduces to the one-mode case: 



1 



n 



max 



E 



max 



{k),n 



X 



Hv^t + yi7) 1 

< TV, + _ 

2 - 2 



(fc),n-r r(/£) ^(k),7l 



(39) 



where we have rewritten the input energy constraint in two 
steps 



W.n ^ y /2 < AT^. + 1/2, 



(40) 



(k) (k) 

and the maximization is over both the CMs V^^. , 

and over the positive integers Nk under the constraint 
EkNk/n^N. 

In conclusion, for quantum memory channels that can be 
unraveled, the calculation of the Holevo function has been re- 
duced to the case of independent, but not identical, one-mode 
channels, each characterized by the matrices x'l'^^'", Y^''^'"' 
and a maximal number Nk of excitations per mode. The only 
ingredient that mixes the one-mode channels is the constraint 
on the total number of excitations. 



Since the Holevo function is invariant under unitary trans- 
formations, the transformation (|32] | preserves the Holevo 
function, i.e. 



1$^ 



(35) 



Moreover, since the encoding Gaussian unitary preserves the 
form of the energy constraint, we have 



tr [E„iVin + V,)El] - tr [Vir, + V,] 



(36) 



which is satisfied if Ej^En = 1, i.e. if the matrix En is orthog- 
onal. Symplectic matrices that are also orthogonal constitute 
a subgroup whose elements are characterized by the form in 
Eq. (O. 

For a memory channel that can be unraveled, it is natural to 
restrict to Gaussian input ensembles such that 



k—l,...n 



and 



EnV,El= Vl 



(fc) 



(37) 



(38) 



A. Examples 

The properties defining the class of Gaussian memory chan- 
nels that can be unraveled are rather peculiar, however several 
relevant models of Gaussian channels belong to this class. In 
this section we review some examples of Gaussian memory 
channels that can (or cannot) be unraveled. 

Lossy Bosonic memory channel. We refer to the general 
model of lossy bosonic Gaussian channel with memory that 
has been introduced in fr?]. Upon 71 uses of the channel, 7i 
input modes are mixed with a corresponding set of n environ- 
mental modes at a beam splitter with transmissivity rj. In the 
Heisenberg picture the canonical field operators transform as 



Qk 
Pk 



V?? Qk + V^-rj Qk 



(41) 
(42) 



/c— l,...n 



where {Qk, Pk}k=i....n are the field operators of the environ- 
mental modes. The environmental modes are in a correlated 
Gaussian state, which is characterized by a CM Vcnv with non- 
vanishing off-diagonal terms coupling different modes. The 
memory channel is associated with the matrices X„ = ^yr^^2n, 
Y,i = (1 — 77) Vcnv Memory effects appear in the channel 
if the n-mode environment is in a non-factorized state. A 
factorized state is characterized by a block-diagonal n-mode 
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CM, i.e. Vcnv = ^^fc- Remarkable examples of non- 

factorized states are the multimode entangled states, which 
belong to the family of multimode squeezed states. To ad- 
dress the problem of memory unraveling, we notice that the 
CM Vcnv can be diagonalized by a 2n x 2n orthogonal ma- 
trix O, i.e. OVcnvO^ = ©JJ^i Vk- Hence one could identify 
Dn := O, En := . Notice that the orthogonal matrix 
preserves the trace, hence it preserves the energy constraint. 
However, a Gaussian unitary is represented with a symplectic 
matrix, hence the given orthogonal matrix represents a physi- 
cal transformation only if it is also symplectic. In conclusion, 
a lossy bosonic memory channel can be unraveled if the CM 
of the environment can be diagonalized by a linear transfor- 
mation which is both symplectic and orthogonal. A character- 
ization of this class of environment CMs is presented in 1 18]: 
pure Gaussian states and squeezed thermal states belong to 
this class. 

Additive noise channel. For this class of channels, the field 
operators transforms according to the Heisenberg picture map: 

qk qk + tk, (43) 

Pk Pk + Uk, (44) 

where tk, Uk are classical stochastic variables. The channel is 
Gaussian if the noise variables are Gaussian distributed, with 
a CM Vc\. The subscript 'cl' refers to the fact that Vc\ is a 
classical CM, i.e. it is only subject to the conditions of being 
symmetric and positive semi-definite. Differently from the 
quantum CM, it needs not obey the Heisenberg uncertainty 
relations as expressed by Eq. (|6]l. Memory effects arise when 
the Vc\ has nonvanishing off-diagonal terms coupling different 
channel uses. The matrices associated to n uses of the channel 
are Xn = ^2n, and Yn = Vc\. The memory channel can hence 
be unraveled if there is a 2n x 2n matrix S, being both sym- 
plectic and orthogonal, such that SnVc\Sj^ — Vk- The 
encoding and decoding Gaussian unitaries which unravel the 
memory channel are hence chosen as £)„ := Sn, En := Sj^. 
Notice that since Sn is orthogonal, — Sn ^ ■ Examples of 
additive noise channels with memory that can be unraveled 
were studied in lfl9ll20ll . The conditions on the CM Vc\ can be 
easily obtained as it is done for the lossy channel in 1 18], with 
the only difference that Vc\ needs not obey the uncertainty re- 
lations. 

Inter-symbol interference channels. In this family of quan- 
tum channels memory effects come from the fact that the sig- 
nals at different channel inputs do interfere at the channel out- 
put, while in the previous examples they are caused by noise 
correlations. For such a case, one can write the Heisenberg 
picture transformations acting on the field operators as 

h 

Pk ^ Y.{<'^'i'^+^'^^''PH + ^^^^'Qk + ntPPn] , 

h 

where the operators {Qfe,Pfc} correspond to environmental 
modes, and the matrices M^??), M(«p), . . . , ^W)^ ^ipP) sat- 
isfy proper conditions |4]. To describe this channel, we work 



in the representation (O. In this representation, assuming that 
the environmental modes are in a Gaussian state character- 
ized by the CM Vcnv, the matrices associated with the memory 
channel are 

^n-y j^(pq) j^(pp) j , (45) 

and 

" ~ 1^ n^pQ) n^pp) ) n^pQ) n^pp) ) ■ ^ ' 

An instance of this kind of model was considered in ||2lll , in 
which M(««) = M^PP\ M^'JP) = M^^') = O, N^^'^) = n(pP)^ 
^{qP) — i^(pQ) = and the environment is in the vacuum 
state, i.e. V^nv = l2n/2. As is shown in |.21], such a channel 
can be unraveled by performing the singular value decompo- 
sition of the matrix X„. 



IV. OPTIMIZATION UNDER SYMMETRIES 

In the cases of both memoryless and memory quantum 
channels, one can pose the question of finding the optimal in- 
put ensemble, i.e. the most robust one under the action of the 
noisy channel. In the memoryless setting, this is related to the 
issue of additivity of the Holevo information: if the Holevo in- 
formation is additive the optimal input ensemble constitutes of 
states which are separable among different channel uses; oth- 
erwise, for channels with a non-additive Holevo information 
the optimal input ensemble is made of entangled states. For 
memory channels it has been observed, in the cases of both 
discrete |2] and continuous |3] variables, that entangled input 
states may be optimal to maximize the Holevo information. 
In models for quantum channels with memory, is customary 
to introduce a memory parameter, used to quantify the mem- 
ory in the channel, which vanishes in the memoryless limit. It 
may happen that the optimal input ensemble constitutes of en- 
tangled states when the memory parameter is above a certain 
threshold. In this case one says that the memory channel ex- 
hibits a transitional behavior. In the case of discrete variables, 
several models exhibit a finite value of the memory threshold 
1 2], while for continuous-variable models the threshold value 
may vanish 1 3, 20], i.e. entangled input states are optimal even 
for arbitrary small, but not zero, values of the memory param- 
eter. 

Restricting to the case of bosonic Gaussian channels that 
can be unraveled, here we introduce a criterion to decide 
whether the Holevo information, restricted on Gaussian input 
ensembles and under input energy constraint, is optimized by 
separable input states. We formulate the following: 

Criterion 1 Given a Gaussian memory channel, represented 
by the sequence (0, X„, Yn), and provided that it can be un- 
raveled, a necessary condition for the optimality of entangled 
input states is the non invariance under phase rotation. 

Proof of Criterion [T] For a Gaussian memory chan- 
nel that can be unraveled, it holds (0,X„,y„) ~ 
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(0= ®k=i,...n X?^'"^ ®k=u...n ^'")- Let us assume, by 
contradiction, that the one- mode channels (0, y/'')'") 
are invariant under phase rotation, i.e. 



where the 2x2 matrix Ri{6) represents a phase rotation 



(47) 
(48) 



Ri{e) := 



cos ( 
sin^ 



— smf 
cos 9 



Matrices that have this symmetry are scalar, i.e. 

(fc),« f xC^)'" 





yC^')'' 



(49) 

(50) 
(51) 



It follows that the solution of the optimization problem 
is also invariant under phase rotation and is given by the ma 
trices: 



1/2 
1/2 



(fc) 



Nk 
Nk 



(52) 



where the optimal values of the parameters {Nk}k=i,...n are 
obtained from the maximization problem 



Yri = max — 

^ {N,} n 



'fl 

E{5 [(x(^-)'"f(A^fc + l/2) + ; 
fc=i 



(fe),n 



1/2 



(2.(fc),n)2(^/2)+yW-"- 1/2]} , 



(53) 



where the maximum is under the constraint X]fe=i Nk/n = 
N. 

Let us notice that the matrix Vi'll in i52\ represents the CM 

of coherent states 101 ■ Then, the matrix ©^^i y^C^) " repre- 
sents the CM of the optimal n-mode input state of the dressed 
memory channel, which includes the encoding unitary trans- 
formation. The actual optimal input state is obtained from it 
by undoing the encoding transformation, i.e. 



opt 



En 



Lk=l 



El 



(54) 



However, since the encoding matrix i?„ is orthogonal (and, of 
course symplectic), it follows that 



V 



opt 



e 

k=l 



1/2 
1/2 



(55) 



i.e. the optimal Gaussian inputs are coherent states, which are 
separable among different channel uses. q 



Remark 1 If the decoding symplectic matrix _D„ is also or- 
thogonal, the conditions Ii47\l . Ii48\l can be equivalently formu- 
lated as follows: 



Rn{0)XnRn{9Y — Xn, 
Rn{e)YnRn{ey = Y^, 



(56) 
(57) 



where the matrix Rn{d) — Ri{0) represents a global 

phase rotation on the n modes. 

The equivalence can be readily proved by working in the 
representation (O. It is a consequence of the fact that the sub- 
group of symplectic and orthogonal matrices [having the form 
as in Eq. ( fT3] )1 commutes with phase rotations, which are rep- 
resented by matrices of the following form: 



Rn 



COS Ulr 

sin Olr, 



— smf 
cos 6] 



(58) 



In other words, if both the encoding and decoding symplec- 
tic matrices are orthogonal, the symmetry under phase rota- 
tion can be checked directly on the matrices X^, Yn. 

Furthermore, being the £"„ elements of Lie groups, 
they reduce to identity when the group's parameters reduce to 
zero. Since the latter would characterize the degree of mem- 
ory, we may argue that the transition can only occur at zero 
value of the memory parameters. This is in contrast to what 
happens in discrete quantum memory channels. 

As a consequence of the Criterion [T] entangled Gaussian 
codewords may be necessary for optimizing the Holevo infor- 
mation ( [39] l only if the rotational invariance is broken. Notice 
however that this is a necessary but not sufficient condition. 
Below we are going to present examples for all possible cases. 



V. EXAMPLES 

Let us first consider the case of an additive noise chan- 
nel. The definition of the channel and its basic properties are 
briefly recalled in Sec. IIII Al Let us recall that for n uses of 
the channel we have Xn = l2n, and Yn = Vc\. Since we 
are considering channels which can be unraveled, we assume 
the existence of an orthogonal and symplectic matrix Sn, such 
that SnVciSl = ©^^1 Vk, with - Sn, En = si. In this 
case, the symmetry condition is verified if the matrix Vc\ is 
symmetric under phase rotation (remark [T]). Two models of 
Gaussian memory channels with Markovian correlated noise 
were studied and characterized in [19, 2Q\. Using the repre- 
sentation ([8]l, the noise CM in ifigll has the following form 



14i 



V O 
O V 



(59) 



which is clearly symmetric under phase rotations ( |58] ). hence 
the optimal Gaussian inputs are separable for this model. On 
the contrary, the noise CM for the model studied in ll20ll has 
the form 



14i = 



V O 

o v 



(60) 
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with V 7^ v. Such a matrix is not symmetric under phase 
rotations ( |58] l and, as shown in [20,1 . the optimal input states, 
when restricted to Gaussian states, are entangled. 

The case of lossy bosonic memory channel |17] is analo- 
gous to the additive channel. Its basic properties are reviewed 
in Se dnTAl In this case X„ = Y^l2„,andF„ (1 — ?/)Knv 
A model of memory channel belonging to this family has been 
studied in [ 1 8J . Using the representation the environmen- 
tal CM in llisll has the form 



1 



(61) 



where M is a symmetric matrix of size n, and T, s are two 
positive parameters. The parameter s quantifies the amount of 
memory in the channel: For s — the environmental state is 
an uncorrected thermal state, while it is entangled for s ^ 0. 
For all s ^ 0, the matrix Yn is not symmetric under phase ro- 
tation and, as shown in ifisll . the optimal Gaussian input states 
are entangled among different channel uses. 

The general case of an inter-symbol interference channel 
is recalled in Sec. IIII Al An example of such a channel was 
studied and characterized in MM, where 



and 



Xri 



y ^If NN^ O 
" ' 2 I O NN^ 



(62) 



(63) 



The channel is invariant under phase rotation, thus, as shown 
in ll2Tll . the optimal input states are separable. To conclude, 
we notice that there are models of Gaussian memory chan- 
nel for which the optimal Gaussian input states are separable, 
even though the channel is not symmetric under phase rota- 
tion. A channel model presenting this feature can be obtained 
by choosing the environment to be in a state of the form 



implying 



and choosing X„ as in ( |62] |. The channel is hence not sym- 
metric under phase rotation, however it is not difficult to show 
that the optimal Gaussian input states are separable. 



VI. CONCLUSION 

In this article we have provided a unified framework for 
some recent results about the performance of quantum Gaus- 
sian memory channels. We have focused in particular on 
the entanglement of optimal input states and we have related 
this issue to the symmetry properties of the channel. More 
specifically we have shown that entangled Gaussian code- 
words might be necessary for optimizing the Holevo informa- 
tion only if the rotational invariance is broken by the channel's 
action. Similar considerations were also done in i fToi l20l I22I1 
for specific channel models. Moreover, for a Gaussian mem- 
ory channel that can be unraveled, we may argue that the tran- 
sition from the optimality of separable states to the optimal- 
ity of entangled states may only occur for vanishing value of 
the memory parameter jlSl 1201 . This is in contrast to what 
happens in discrete quantum memory channels 12(]. However, 
while there investigations have only involved very few chan- 
nel uses, here the analysis has been carried out for arbitrary 
number of channel uses by resorting to a mathematical ma- 
chinery called memory unraveling. That allowed us to trace 
the Gaussian memory channel back to a memoryless one. Sev- 
eral examples have been discussed concerning memory unrav- 
eling as well as transitional behavior. 

We think that the presented results shed light on the mech- 
anisms and the structure of the correlations that lead to an en- 
hancement of the channel performance with entanglement, al- 
though a complete characterization of memory channels tran- 
sition features is still far away. 



(64) 

Acknowledgments 

This work has been supported by the European Commis- 
sion under the FET-Open grant agreement CORNER, number 
(65) FP7-ICT-213681. 



[1] R. G. Gallager, Information Theory and Reliable Communica- 
tion (Wiley, New York, 1968). 

[2] C. Macchiavello and G. M. Palma, Phys. Rev. A 65, 050301(R) 
(2002); E. Karpov, D. Daems and N.J. Cerf, Phys. Rev. A 74, 
032320 (2006); D. Daems, Phys. Rev. A 76, 012310 (2007); 
F. Caruso, V. Giovannetti, C. Macchiavello and M.B. Ruskai, 
Phys. Rev. A 77, 052323 (2008). 

[3] G. Ruggeri, G. Soliani, V. Giovannetti and S. Mancini, Euro- 
phys. Lett. 70, 719 (2005); G. Ruggeri and S. Mancini, Quan- 
tum Inf. Comput. 7, 265 (2007); O. Pilyavets, V. Zborovskii and 



S. Mancini, Physs Rev. A 77, 052324 (2008). 
[4] A. Ferraro, S. Olivares and M. G. A. Paris, Gaussian states 

in quantum information (Bibliopolis, Napoli, 2005); eprint 

arXiv:0503237 [quant-ph]; S. L. Braunstein, P. van Loock, Rev. 

Mod. Phys. 77, 513 (2005). 
[5] R. Simon, N. Mukunda, and B. Dutta, Phys. Rev. A 49, 1567 

(1994). 

[6] A. S. Holevo, R. R Werner, Phys. Rev. A 63, 032312 (2001). 

[7] M. M. Wolf, Phys. Rev. Lett. 100, 070505 (2008). 

[8] F. Caruso, J. Eisert, V. Giovannetti and A. S. Holevo, New J. 



9 



Phys. 10, 083030 (2008). 
[9] A. S. Holevo, arXiv:quant-ph/060705r 

[10] A. S. Holevo, IEEE Trans. Inf. Theory 44, 269 (1998); B. Schu- 
macher, M. D. Westmoreland, Phys. Rev. A 56, 131 (1997). 

[11] V. Giovannetti, S. Guha, S. Lloyd, L. Maccone, J. H. Shapiro, 
H. P Yuen, Phys. Rev. Lett. 92, 027902 (2004). 

[12] M. B. Hastings, Nature Physics 5, 255 (2009). 

[13] G. Bowen, S. Mancini, Phys. Rev. A 69, 012306 (2004). 

[14] D. Kretschmann, R. R Werner, Phys. Rev. A 72, 062323 (2005). 

[15] V. Giovannetti, S. Guha, S. Lloyd, L. Maccone, J. H. Shapiro, 
Phys. Rev. A 70, 032315 (2004). 

[16] A. S. Holevo, R. R Werner Phys. Rev. A 63, 032312 (2001); A. 
S. Holevo, M. Sohma, O. Hirota, Phys. Rev. A 59, 1820 (1999). 



[17] V. Giovannetti and S. Mancini, Phys. Rev. A 71, 062304 (2005). 
[18] C. Lupo, O. V. Pilyavets, S. Mancini, New J. Phys. 11, 063023 
(2009). 

[19] C. Lupo, L. Memarzadeh, S. Mancini, Phys. Rev. A 80, 042328 
(2009). 

[20] J. Schafer, D. Daems, E. Karpov, N. J. Cerf, Phys. Rev. A 80, 
062313 (2009). 

[21] C. Lupo, V. Giovannetti, S. Mancini, Phys. Rev. Lett. 104, 
030501 (2010). 

[22] N. J. Cerf, J. Clavareau, C. Macchiavello, J. Roland, Phys. Rev. 
A 72, 042330 (2005). 



